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The diffusion Monte Carlo method is applied to describe 
a trapped atomic Bose-Einstein condensate at zero tempera- 
ture, fully quantum mechanically and nonperturbatively. For 
low densities, n(0)a^ < 2 • 10~^ ['^■(0): peak density, a: s- 
wave scattering length], our calculations confirm that the 
exact ground state energy for a sum of two-body interac- 
tions depends on only the atomic physics parameter a, and 
no other details of the two-body model potential. Correc- 
tions to the mean-field Gross-Pitaevskii energy range from 
being essentially negligible to about 20 % for TV = 2 — 50 
particles in the trap with positive s-wave scattering length 
a — 100 — 10000 a.u.. Our numerical calculations confirm 
that inclusion of an additional effective potential term in the 
mean-field equation, which accounts for quantum fluctuations 
[see e.g. E. Braaten and A. Nieto, Phys. Rev. B 56, 14745 
(1997)], leads to a greatly improved description of trapped 
Bose gases. 
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I. INTRODUCTION 

Since the first achievement of Bose-Einstein condensa- 
tion in trapped atomic vapors in 1995 O, these systems 
have received increased attention from both experimen- 
tal and theoretical efforts. Theoretical studies of these 
inhomogeneous gases are primarily based on the mean- 
field Gross-Pitaevskii (GP) equation for the condensate 
wave function ||]. This equation can also be viewed as 
a variant of the Hartree-Fock (HF) approximation [p|. 
It is thus fundamentally important for our understand- 
ing of this many-boson system to ascertain the validity 
of the mean-field description and the importance of par- 
ticle correlations. In principle, these questions can be 
investigated through the use of a complete many-body 
basis instead of making a single particle approximation. 
Some studies have in fact considered the T — con- 
densate at the level of the random phase approximation 
(RPA), with a few going still further to treat the particle- 
particle correlations at the level of configuration interac- 
tion (CI) Pfl. For many particles, however, a full CI 
calculation exceeds current computational capabilities. 
Consequently we adopt the diffusion Monte Carlo (DMC) 
method for the calculations reported in this paper. 

In the mean-field description, the particle interactions 
enter solely through an effective mean-field potential, 
which is proportional to the s-wave scattering length 



a [Ep-nS]. In contrast, a full description treats a many- 
body potential surface and allows nonseparability of the 
many-body wave function. The first key question of our 
study is therefore: Can the s-wave scattering length ap- 
proximation properly describe an inhomogeneous Bose 
gas, or does the actual form of the two-body inter- 
action potential necessarily come into play? To test 
the validity of the shape independent atom-atom po- 
tential (frequently expressed as a (5-function potential 
with s-wave scattering length a [M), we solve the TV- 
body Schrodinger equation for an inhomogeneous Bose 
gas for various different two-body potentials that gen- 
erate identical a using the DMC method for densities 
n{0)a^ < 2 ■ 10~^. We find that different potentials 
produce indistinguishable total ground state energies E, 
indicating that the lowest many-body Schrodinger en- 
ergy eigenvalue is independent of the shape of the two- 
body potential. This result is in a greemen t with predic- 
tions based on an expansion in y^n{0)a^ for low densi- 
ties 0,0. 

However, we do observe differences between the accu- 
rate many-body ground state energy and that obtained 
by solving the GP equation. A second motivation of our 
study is therefore to assess the validity of the mean-field 
description as a function of a and iV, and as a function of 
the peak density n(0)a^. Furthermore, we test the valid- 
ity of a modified GP equation that accounts for quantum 
fiuctuations [T^jl^. In addition to the ground state en- 
ergy, we also show that the condensate density can be 
significantly affected by correlations. These calculations 
support our conclusion that the modified GP equation 
[see Eq. (S)] greatly improves upon the commonly used 
mean- field treatment. 

Theoretical studies of correlation effects in trapped 
condensates have already been proposed g,^,0-g7|. 
These approaches include approximate zero tempera- 
ture variational Monte Carlo studies |2^, essentially ex- 
act finite temperature path integral Monte Carlo stud- 
ies |5|,|6|, po|,pi[ as well as correlated basis function ap- 
proaches [^2[ , and perturbative schemes |p7|- p9| , p7| . Here, 
we use essentially exact diffusion Monte Carlo techniques 
to directly calculate correlation effects at zero tempera- 
ture. 

Section || introduces the many-body Schrodinger equa- 
tion of a trapped Bose gas and the numerical treatment 
applied to solve this equation. Results and their inter- 
pretation are presented in Section [II. Section IV sum- 
marizes this paper. 



II. SYSTEM AND NUMERICAL TECHNIQUES 

The many-body Schrodinger equation for a condensate 
of A'^ mass m bosons in a spherical trap, centered at the 
origin with trapping frequency LUho, is given by 
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tjj = Eip. 
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Here V{nj) = V{\r 



j\) denotes the two-body inter- 



action potential, and fi = (x^, y^, Zi) is the Cartesian po- 
sition vector of atom i relative to the trap center, m is 
taken to be m(^'^Rb), and ujho to be 2tt x 77.78 Hz. In 
the following we use harmonic oscillator units for ener- 
gies {hiVho) and length {aho = [h/{mujho)]^^'^)- The en- 
ergy per particle for the non-interacting case, V{r) = 0, is 
then E/N = 1.5 hujho, and, for example, the *^Rb s-wave 
triplet scattering length (a « 100 a.u.), a = 0.00433 a^o 
at that frequency. 

The DMC method [E8[ (and references therein), here 
implemented with importance sampling and a descen- 
dant weighting scheme, derives a solution to the many- 
body ground state Schrodinger equation, Eq. (||), for a 
given model potential surface. This solution is essen- 
tially exact, apart from statistical uncertainties. The re- 
sulting ground state energy shall be denoted by Edmc 
in the following. Ideally, one would solve Eq. p) for 
the best known ®^Rb-*^Rb two-body potential V{r) (e.g. 
Rcf. [^). However, the large number of bound states 
in this potential makes its use in a many-body calcula- 
tion extremely challenging, if not infeasible (see below). 
Instead, we have carried out tests using three different 
two-body model potentials A-C to test how the energet- 
ics depend on the actual form of the two-body potential: 

A) a purely repulsive hard core potential with hard core 
radius a, V{r) = oo for r < a and V{r) = for r > a; 

B) another purely repulsive potential with parameters 
d > and tq, namely V{r) — dcosh~^(r/ro); and C) a 
sum of an attractive and a repulsive Gaussian, leading 
to a potential that exhibits a minimum with negative en- 
ergy. The potential parameters are adjusted such that i) 
a = 0.00433 aho (potentials A-C), and ii) a — 0.433 a^o 
(potentials A and B). 

Potential C requires a few more remarks. In contrast 
to potentials A and B, which are purely repulsive, poten- 
tial C is repulsive at short range and attractive at long 
range. The usage of this potential in a DMC calcula- 
tion could therefore potentially lead to a lowest energy 
state that describes a cluster rather than a condensed 
state. However, the two-body potential is such that it 
does not support a two-body bound state; nevertheless, 
many-body bound states might exist. If this two-body 
potential supports a many-body bound state, then our 
simulation should converge to a cluster state rather than 
to the metastable condensed state at infinite imaginary 



simulation time. At finite time, however, we find empiri- 
cally for positive a that our DMC simulation samples the 
condensed and not the cluster state. 

To test our DMC code we carry out a separate exact 
solution to Eq. (|l]) for two particles, iV = 2, by sep- 
arating out the center of mass motion and solving the 
one-dimensional radial Schrodinger equation using stan- 
dard numerical techniques. For potential A we find, for 
example, E = 3.00345 huiho for a — 0.00433 aho, and 
E = 3.3827 hcuho for a = 0.433 aho in excellent agree- 
ment with our DMC results, Edmc — 3.00346(1) hcuho 
and Edmc — 3.3831(7) hujho, respectively. (The statis- 
tical error of the DMC energies is given in parenthesis.) 
Since the DMC algorithm generalizes straightforwardly 
as a function of iV, the above agreement is a convincing 
test of our code. Furthermore, we carefully checked the 
independence of our results on the exact shape of the trial 
wave function and on the time step used for the propaga- 
tion in imaginary time. The next section presents results 
from our DMC calculation and compares them with GP 
theory. 



III. RESULTS AND INTERPRETATION 

Columns 2-4 of Table | show the DMC energies Edmc 
for iV = 3 -- 20 for the three two-body potentials A, B 
and C with i) a = 0.00433 aho (top) and ii) a = 0.433 aho 
(bottom). For the parameter range considered in Table |, 
which corresponds to n{Q)a^ < 2 • 10"'^ (see below), the 
ground state energy Edmc of the inhomogeneous gas for 
different V{r) with identical two-body scattering length 
a is independent of the shape of V{r); and thus Edmc 
depends only on a to within our statistical uncertain- 
ties. For the homogeneous Bose gas, Giorgini et al. pC| ] 
find a small dependence of the ground state energy on 
the shape of the two-body potential in the intermediate 
density regime, na^ « 10~^ {n: number density of ho- 
mogeneous gas); however, their energy depends only on 
a in the low density regime. Recent extensive studies of 
trapped two- and three-particle condensates as a func- 
tion of the scattering length a and the trap frequency 
Uho reveal dependences of the ground and excited state 
energies on the exact shape of the two-body potential in 
the high density limit. However, these studies will be 
presented elsewhere. 

To ascertain the validity of the widely used mean-field 
approach, we now compare our many-body DMC ener- 
gies with those obtained by solving the GP equation. 
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Here, <I>Gp(r) denotes the ground state orbital (normal- 
ized to 1), and ecp the orbital energy; the total energy 
Eqp can be obtained through evaluation of the energy 



functional Ecp[^gp] 0- Notice that ecp can also be as- 
sociated with the chemical potential of the system. Both 
ecp and Egp depend only on {N — l)a, rather than on 
a and N separately. The formulation of the GP equation 
in terms of [N — l)a rather than Na follows from number 
conserving Schrodinger quantum mechanics, using a HF 
initial state js). To compare our DMC energies Eoaic 
with the GP energies Egp, our many-body DMC calcu- 
lations reported below have been performed for potential 
A only, which is the purely repulsive hard core potential. 
Since Edmc is independent of the shape of V{r) (see Ta- 
ble ffl), this justifies our restriction to only one two-body 
potential. Note, however, that our calculations reported 
in the following treat densities 71(0)0"^ as large as « 1. In 
this high density regime we expect the many-body en- 
ergy to depend on the detailed shape of the two-body 
potential. 

Figure |l|(a) compares the ground state energies Egp 
and Ejjmc for A^ = 2 — 50 as a function of the in- 
teraction parameter {N — l)a/aho [Eq. (||)]. Note the 
logarithmic scale for (N — l)a/afio', diamonds and tri- 
angles show our calculated data points, while solid and 
dotted lines connect those data points with identical 
N but different a. The full quantum energies Edmc 
are always equal to or larger than the mean-field en- 
ergies Egp, with [Edmc ~ Egp) /Edmc attaining 20 
to 50 % for the largest scattering lengths a consid- 
ered here, a = 3.9,3.0,2.6,1.7,0.87,0.43 Uho for N = 
2, 3, 5, 10, 20, 50. Figure |l](a) shows the separate depen- 
dences of Edmc on a and N; e.g., for {N — l)a = 10 Uho, 
(Edmc - Egp)/Edmc ~ 0.4 for A = 5 , but instead 
« 0.07 for A^ = 50. Column 5 of Table | summarizes 
Egp for a few selected A^ and a parameters. In conclu- 
sion, the GP equation, Eq. (||), describes the condensate 
accurately for small a. For larger a, however, significant 
departures occur from the DMC results. Later in this 
section we discuss these departures further. 

We refer to the energy difference Edmc ~ Egp as 
a "correlation energy" Ecorr- Note that the definition 
of Ecorr introduced here differs from that often used in 
quantum chemistry, where Ecorr is defined as the en- 
ergy difference between the CI energy Eci and the HF 
energy Ehf for the same two-body potential V{r). In 
our treatment Edmc results from solving the many-body 
Schrodinger equation, Eq. dl|), for the exact wave func- 
tion (= complete basis) with a shape dependent two-body 
potential V{r), whereas Egp results from the mean-field 
equation, Eq. ^), or equivalently, from Eq. (|l|) for a 
shape independent (5-function potential using a HE ini- 
tial state. Thus, Ecorr not only contains basis set ef- 
fects, but also effects due to the reduction of the shape 
dependent interaction potential to a shape independent 
potential that generates an identical scattering length 
a. Note that Edmc > Egp, and therefore Ecorr > 0, 
for the parameter range considered here. The "exact" 
many-body treatment of the particle interaction (using a 
model two-body potential) leads evidently to an increase 
of the mean-field energy Egp , which is somewhat coun- 



terintuitive. Recent calculations for potentials of type B 
with negative potential depth d show, however, that the 
many-body energy can be lower than the GP energy for 
extremely tight traps, i^ho ~ 10 MHz (see, e.g., Fig. 3 
of B ) . In the following we restrict ourselves to the hard 
core potential A. 

Next, we assess a modified GP equation, which has 
been discussed in the literature [^ 18|. Braaten and Ni- 
eto [0 included the effects of quantum fluctuations for 
positive a to generate corrections to the mean-field equa- 
tion. Qualitatively, quantum fluctuations can be viewed 
as representing virtual double excitations of particles out 
of the condensate, which is also related to the ground 
state depletion. Within the Thomas-Fermi (TF) approx- 
imation, these "quantum corrections" produce an addi- 
tional local effective potential term in the GP equation 



that is linear in the assumed small parameter {na ) 
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The total energy EGP.mod relevant to this modifled GP 
equation can be obtained from the energy functional 
EGP,mod\^GP,mo(A- Equation (§) accounts for two-body 
physics through the scattering length a solely; three-body 
effects are not included. In contrast to Eq. (||), which 
only depends on the interaction parameter (A^ — l)a, 
Eq. (H) exhibits separate dependences on A — 1 and a. 
Note that inclusion of second order terms in the TF ex- 
pansion leads to additional non-local terms in Eq. (|^), 
which account for edge effects [|7|. Our estimates sug- 
gest that these non-local terms are small, and we have 
therefore neglected them. 

The term accounting for quantum fluctuations for 
the inhomogeneous gas, the second term in the sec- 
ond line of Eq. (g), is identical to the Huang-Lee- 
Yang correction term for the homogeneous gas [|ll|-{l3[. 
Huang, Lee and Yang showed in 1957 that the many- 
body Schrodinger equation for a Bose gas interacting 
via hard sphere potentials can be equivalently formu- 
lated by replacing the hard sphere boundary conditions 
by a pseudo-potential. This pseudo-potential can then be 
treated in perturbation theory, leading to a lowest order 
ground state energy per particle Eq/N « 27r?i {a/m)n. 
To next order in {na^Y''^, the first order perturba- 
tion energy can be evaluated by summing over momen- 
tum states with wave vector k ^ [|l3| , resulting in 
Eq/N « 2TTh^{a/m)n[l + r^{na^Y''^] for the homoge- 
neous gas. The summation over k states depends cru- 
cially on use of the exact pseudo-potential expansion: 
Use of 4:nh'^{a/ni)S{r) rather than 4:7rh^{a/m)S{r)-^r 
leads to a divergent ground state energy expression. The 
additional term arising in first order perturbation theory 
for the homogeneous gas is identical to the term in the 



inhomogeneous gas accounting for quantum fluctuations 
(using quantum mean-field theory language) |17| , and di- 
rectly shows the contribution of fc ^ (excited) states to 
the ground state energy. Braaten and Nieto ||l^ stress 
that the most important excitations are to orbitals in the 
energy range Eexatation ~ %Tih^{a/m)n. 

The derivation of Eq. (^ is based on two assump- 
tions |1^]: i) the expansion parameter (na'^)^/^ is small, 
and a) 1/ y/na is small compared to the condensate ra- 
dius. To assess the accuracy of these assumptions explic- 
itly, Fig. |l|(b) compares Ecp^mod obtained from Eq. (||) 
with our DMC energy Edmc ■ Consider the N — 2 curve 
first. EjjMC — Eop^rnod IS positivc for small a, switches 
sign at a w 0.2 aho^ reaches its minimum for a ~ 0.87 aho^ 
and becomes positive for a > 1.7 ato- The behavior for 
larger N is similar. The maximal negative deviation for 
the parameter range shown in Fig. ||(b) is about 3 % for 
TV = 5 and a « 0.87 aho- For A^ = 20 and iV = 50, the 
agreement between E^mc ^^nd Eap^od is better than 
3 % for a as large as 0.87 aho (« 20000 a.u.) for N = 20, 
and 0.43 a^o (« 10000 a.u.) for TV = 50, respectively. 
Table ffl summarizes Eqp„ioj^ for a few selected values of 
N and a. 

The modified GP equation, Eq. (^), overcorrects the 
correlation effects for some a for all N considered. Over- 
all, the additional effective potential term in Eq. (y|) im- 
proves the GP treatment substantially for N = 2 — 50. As 
the size of the condensate increases, we expect the next 
higher order correction terms, which are non-local p7[ , 
to become even smaller, and the local correction term 
considered in Eq. (0) to be the dominant one. It should 
be noted, however, that Eq. (H) is derived in the large 
TV regime, where the number of particles contributing to 
quantum fluctuations is small compared to those occu- 
pying the "single particle ground state" |l^. Thus the 
excellent agreement between EGP,mod and Edmc in the 
small N limit might be somewhat fortuitous. 

Figure shows the same data as in Fig. Ma), how- 
ever, now as a function of [N — \){a/ ahoY rather than 
the interaction parameter [N — l)a/aho- AH data points, 
TV = 2—50, lie on a "universal curve" . As we shall discuss 
later, this universal behavior suggests that our calcula- 
tions for small TV and their comparison with mean-field 
and quantum mean-field treatment may directly transfer 
to inhomogeneous gases with larger TV. 

To summarize our ground state energy studies, note 
the different scales of the vertical axes in Figs. ||(a) and 
(b) . A comparison of (a) and (b) shows that a significant 
improvement over GP theory is accomplished through the 
inclusion of quantum fiuctuations [Eq. (y)], as evidenced 
by the better agreement between Ecp^mod and Ejjmc 
than between Eqp and Ep,MC- 

Besides the ground state energy, we also calculate 
structural properties using Eqs. (|l|)-(||). The peak den- 
sity n(0), for example, determines many properties such 
as the lifetime of a trapped condensate, and its exact de- 
termination is therefore of great importance. Figure 
summarizes our results for the radial density per parti- 



cle n{r)/N for TV = 3 (a) and A^ = 10 (b) for two dif- 
ferent scattering lengths a, where J n{r)(Pr — TV. Our 
DMC density per particle nDMc{f)/N [Eq. (|l|)] is shown 
as triangles for a — 0.0433 aho and as diamonds for 
a = 0.433 aho- Note that these data points have sta- 
tistical noise. Dashed lines indicate the GP density per 
particle ncp{r)/N [Eq. (0)], while dotted lines show the 
modified GP density per particle ncp^modir)/N [Eq. (||)]. 
For the smaller scattering length, a — 0.0433 aho, 
all three densities per particle, nDMc(j')/N, nQp{r)/N 
and ncp,mod{f)/N, agree well. In fact, for TV = 3 
the modified GP density and the GP density are in- 
distinguishable. For the larger scattering length, a — 
0.433 aho, ncp,mod{'i')/N coincides with the DMC density 
ri£)Afc'(r)/TV, whereas the GP equation overestimates the 
peak density by about 20 % (TV = 3), and 35 % (TV = 10), 
respectively. For comparison, the GP treatment underes- 
timates the energy for this scattering length by only 4 % 
(10 %), whereas the modified GP equation overestimates 
the energy by only about 2 % (2 %) for TV = 3 (TV = 10). 
The calculation of the DMC density udmc {f) , and es- 
pecially the peak density nuMci^), is niore time inten- 
sive than that of the DMC energy Eomc ■ To replot our 
data from Fig. we therefore approximate noMc{^) by 
nGP,mod{^), which is easy to calculate. Figure ^ shows 
the fractional energy correction {Edmc ~ Eqp)/Eumc 
on a log-log scale as a function of 7iGP,mod(0)a'^ for 
lO^'* < nap,mod{^)o?' < 1 (sohd and dotted lines). In 
addition, the dashed line shows the fractional correction 
to the GP energy predicted within quantum mean-field 
theory using the TF (large TV) limit, x/(l + x), where 



X = 7/8y^Tm{0)a^ iO-Ill- The larger N, the smaller 
the difference between the fractional correction to the GP 
energy based on our DMC calculation for small TV and the 
asymptotic large TV limit (dashed curve). Figure seems 
to suggest that the predicted behavior based on quantum 
mean-field theory within the TF limit [p yi7| , p^ is indeed 
correct, and gives an accurate estimate of the error of 
the GP equation as a function of the density n{0)a'^ for 
large TV. Our calculations may therefore be viewed as a 
first explicit numerical test of the application of quantum 
mean-field theory to trapped Bose gases. 



IV. SUMMARY 

This paper presents full quantum calculations for an 
inhomogeneous Bose gas that serve as a stringent test of 
mean-field theory. Our calculations suggest that conden- 
sates are well described through only the s-wave scat- 
tering length a. Other characteristics of the two-body 
potential such as the effective range parameter seem to 
be neglegible for the a and TV parameter range con- 
sidered here. These findings for the inhomogeneous 
gas are in agreement with quantum mean-field predic- 
tions [ p7| , p^p8| , and recent DMC studies for the homo- 
geneous Bose gas [pffl . Our studies also suggest that the 



first order correction term [Eq. (|^)] leads to a great im- 
provement upon the GP equation, [Eq. (g)], for large a. 
For a spherical trap with, for example. A'' = 5000 par- 
ticles and a = 0.0433 a^o (1000 a.u.) the GP and the 
modified GP energies differ noticeably, at the 5.9 % level: 
Egp/N = 9.239 hiVho, and Eap^mod/N = 9.786 hcvho- 
The densities for this condensate differ at the 15 % level: 
nGp{0)a^ = 1.9-10-3 and nGP,mod{0)a^ = 1.7- lO'^. We 
expect the excitation frequencies to change by an amount 
comparable to that of the energy [pip^ ] . Very recently, 
stable condensates with large scattering length a — 
10000 a.u., corresponding to a density n{0)a^ — 10"^, 
have been obtained experimentally in a non-spherical 
trap |3^ . 

The studies presented here are performed for a two- 
body potential with positive scattering length a, which 
supports no two-body bound state. It will be interest- 
ing to extend these studies to two-body potentials with 
negative a. Furthermore, the interplay between recom- 
bined molecular (i.e. "snow flake") and BEC-type states, 
completely neglected in this work, promises rich physics. 
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FIG. 1. Energy difference Edmc — Egp normalized by 
Edmc (a) and Edmc - Eop.mod normalized by Edmc (b), 
respectively, as a function of the GP interaction parameter 
(A'^ — l)a/aho- Diamonds and triangles show data points 
resulting from solving Eqs. ([l|)-te|), and solid and dotted 
lines connect data points for identical TV to guide the eye, 
N = 2 — 50. Note the logarithmic (A'^ — l)a/aho scale. 



FIG. 2. Energy difference Edmc — Egp normalized by 
Edmc, also shown in Fig. nj(a), however, now as a func- 
tion of (A'' — l){a/aho)^ rather than the interaction param- 
eter (A*' — l)a/aho- AH data points, N = 2 — 50, line up on 
a "universal curve". Note the logarithmic (A'^ — \){a/aho)^ 
scale. 
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FIG. 3. DMC density per particle noMc{i')/N togetiier 
with GP density per particle nGp{r)/N and the modified 
GP density per particle nGP,mod{r) / N iov N = 2, (a) and 
A'^ = 10 (b) for two different scattering lengths a. Triangles 
show the DMC density for a — 0.0433 aho, diamonds that for 
a = 0.433 aho- Dashed lines show the GP and dotted lines the 
modified GP densities per particle. The statistical error of the 
DMC data is smaller than the size of the symbols at large r, 
and about twice the size of the symbols near the trap center 
(small r). Note nGp{r)/N and nGp,mod{r)/N are indistin- 
guishable on the scale shown for A'^ = 3 and a = 0.0433 aho- 



TABLE I. Ground state energies Eomc, columns 2-4, in 
units of huiho for three different two-body potentials A-C for 
A'' = 3 — 20 particles in the trap; top for a = 0.00433 aho 
(potentials A-C), and bottom for a = 0.433 aho (potentials 
A-B). The statistical uncertainty is given in parenthesis. For 
comparison columns 5 and 6 contain the GP energy Eqp and 
the modified GP energy EGP,mod, respectively. 
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FIG. 4. Energy difference Edmc — Eqp normalized by 
Eomc, also shown in Fig. nl(a) and H, however, now as a 
function oi nGp,mod{0)a^ (solid and dotted lines) for densities 
10~* < nGP,mod{0)a^ < 1 together with the fractional correc- 
tion to the GP energy calculated within quantum mean-field 
theory using the TF approximation (dashed line, see text). 
Note the logarithmic scale of both axes. 



N 


potential A 


potential B 


potential C 


Egp 


EGP,mod 


3 


4.51036(2) 


4.51037(2) 


4.51035(4) 


4.51032 


4.51032 


5 


7.53443(4) 


7.53439(6) 


7.53441(10) 


7.53432 


7.53434 


10 


15.1537(2) 


15.1539(2) 


15.1536(8) 


15.1534 


15.1535 


20 


30.640(1) 


30.639(1) 


30.638(2) 


30.638 


30.639 


3 


5.553(3) 


5.552(2) 




5.329 


5.611 


5 


10.577(2) 


10.574(4) 




9.901 


10.772 


10 


26.22(8) 


26.20(8) 




23.61 


26.84 


20 


66.9(4) 


66.9(1) 




57.9 


68.5 



